Nucleation pathway of core-shell composite nucleus in size and composition space and 

in component space 
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The kinetics of nucleation of a core-shell composite nucleus that consists of a core of stable final 
phase surrounded by a wetting layer of intermediate metastable phase is studied using the kinetic 
theory of binary nucleation not only in the size and composition space but also in the component 
space. The steady-state solution of the Fokker-Planck equation is considered. Various formulas 
for the critical nucleus at the saddle point as well as for the post-critical nucleus are derived. 
The kinetics of nucleation at the saddle point is more appropriately characterized in the size and 
composition space, while the kinetics of the post-critical nucleus is more appropriately described 
in the component space. Although both the free-energy landscape and the reaction rates play 
decisive role to determine the kinetics of nucleation at the saddle point, the details of the free 
energy landscape are irrelevant to the kinetics of the post critical nucleus. 

PACS numbers: 64.60.Q- 
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I. INTRODUCTION 

Nucleation is a very basic phenomena which plays a vi- 
tal role in various material processing applications rang- 
ing from steel production to food and beverage indus- 
tries Recently, researchers have focused on the nu- 
cleation of complex materials [l| . The nucleation of such 
complex materials can also be complex and often involves 
intermediate metastable phases [2|-[6| , which appears dur- 
ing the course of nucleation and growth as predicted from 
the Ostwald's step rule 0). Then, the critical nucleus 
often has a core-shell structure with stable final phase 
surrounded by an intermediate metastable phase. 

Such a core-shell nucleus appears in various circum- 
stances. It is well known, for example, the nucleation 
of protein crystal proceeds through the core-shell type 
nucleus with the final stable crystal surrounded by the 
metastable dense solution Q. The model calculation Q 
using the Monte Carlo simulation with a simplified inter- 
molecular interaction revealed that the critical nucleus 
at the saddle point corresponds indeed to the core-shell 
structure. Subsequent numerical simulation using the 
Lennard- Jones system [8| and a model calculation us- 
ing the capillarity approximation 9] confirmed the core- 
shell structure of critical nucleus which corresponds to 
the saddle point of the free-energy landscape. A similar 
composite nucleus is predicted even for a simple metal 
like aluminum JlOJ using molecular dynamics simulation. 
Experimental evidence of such a core-shell structure of 
critical nucleus is observed not only in protein crystal- 
lization but also in colloidal crystallization [lj [l2j • 

A similar composite nucleus with core-shell structure 
is considered in the problem of delquescence [U, E3], 
where the condensation of liquid from supersaturated va- 
por occurs on a soluble core. This core-shell structure 
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also appears in various problems such as semiconduc- 
tor nano-crystals [IH, polymer crystallizations and 
nano-clusters of alloys |17| . It is also used as a model 
of linked-flux nucleation or partitioning transformation 
when the interface-limited growth and diffusion of ma- 
terial is coupled [l8r]2lj|. In this model, a core of the 
final stable phase is surrounded by a shell of the dense 
pre- nucleus environment |18| . 

In the previous paper [22|, we have studied the nucle- 
ation flux of composite nucleus with core-shell structure 
at the saddle point in the free-energysurface using the 
theory of decay of metastable phase |23l - l2q and of mul- 
ticomponent nucleation [27l l28j. In this paper, we will 
supplement our previous study (22[, and will study the 
nucleation pathway for the critical as well as the post- 
critical nucleus with core-shell structure. We first recap- 
ture the Fokker-Planck or the Zeldovich-Frenkel equa- 
tion 0, [2t| for the composite nucleus [13] from the Mas- 
ter equation (Section [TTJ) by regarding the composite nu- 
cleus as a fully phase-separated two-component binary 
system. Then the nucleation flux will be characterized 
more appropriately using the size-composition represen- 
tation instead of the two-components representation (22j 
(Section ing . The nucleation pathway of the post-critical 
composite nucleus will also be studied using the Zeldovich 
relation [3(| Finally, section HVl will contain the con- 
clusion of the study, indicating the outcome as well as 
suggestions on further research in the field. 



II. FOKKER-PLANCK EQUATION FOR THE 
COMPOSITE NUCLEUS 



In order to study the nucleation kinetics of the com- 
posite nucleus with core-shell structure, the model shown 
in Fig. Q] has been used [9j, l22| . The model consists of a 
core of the stable new phase (number of molecules ni) 
surrounded by an intermediate metastable phase (num- 
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ber of molecules 712) nucleated in the metastable parent 
phase. Kashchiev and coworkers [32],[33| have also consid- 
ered this model as a model of nucleus when there exists 
an intermediate metastable phase. 

This core-shell nucleus is considered to form by the 
two-step mechanism: First, the nucleus of the metastable 
intermediate phase appears within the metastable par- 
ent phase. Next, the core of the stable new phase starts 
to nucleate inside the nucleus of the metastable inter- 
mediate phase to form the core-shell structure shown in 
Fig- HI However, in contrast to the naive expectation_|33[ 
the computer simulation Q , mean- field calculation Qas 
well as the experimental result in colloidal crystalliza- 
tion [l2[ have shown that this two-step mechanism is in 
fact the single-step nucleation characterized by a single 
activation energy and, therefore, by a single saddle point. 
This saddlepoint corresponds to the core-shell critical 
nucleus @, H El of the specific size and composition. 
Therefore, we will consider the simplest example of the 
core-shell nucleus of a single-component system. More 
complex scenarios will be expected for multi-component 
system such as the problem of deliquiescence (13 . Il4j and 
of semiconductor nano-crystals [l5| . 
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FIG. 1. A core-shell critical nucleus model that consists of 
a stable new phase (number of molecules m) surrounded by 
an intermediate metastable phase (number of molecules 722) 
nucleated in the metastable parent phase [t| ■ Transformation 
rates k + and k~ are the reaction rates between the stable new 
phase and the metastable intermediate phase. Attachment 
rates a + and a~ are the reaction rate from the metastable 
parent phase to the intermediate metastable phase. 

Transformation rates n + and k~ in Fig. Q] are the 
reaction rates between the stable new phase and the 
metastable intermediate phase. Attachment rates a + 
and a~ are the reaction rate from the metastable parent 
phase to the intermediate metastable phase. The forma- 
tion of stable final phase (core) is assume to occur only 
through the transformation of surrounding metastable 
phase (shell). 

The Master equation for the time-dependence of the 
concentration of clusters f (ni,n 2 ,t) that consists of n\ 
molecules of the stable phase in the core and n 2 molecules 
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of the intermediate metastable pha se in the surrounding 
shell is written generally [H, l29l. l34j in the form 

df(n u n 2 ,t) + 

= or {ni,n 2 - 1) / [ni,n 2 - l,t) 

- [a + (m, n 2 ) + a~ n 2 )] / (ni, n 2 , t) 
+ a~ (nx,n 2 + l)/(ni,ra 2 + l,t) 
+ k + (ni - 1, n 2 + 1) / (ni - I, n 2 + 1, t) 

- [k + (711,712) + K_ (^1,712)] / (ni,n 2 ,t) 
+ K~ (ni + 1, 712 — 1) N (ni + 1, 712 — 1) , 

(1) 



Using the detailed balance condition 
k~ (ni + I, n 2 — 1) = k + (711, n 2 ) 



f eq (ni,n 2 ) 



/eq ("1 + l,7l 2 - l)' 



- / ,11 + / \ /eg (ni,n 2 ) 
a [ni, n 2 + I) = (m, n 2 ) 



/ eq (ni,n 2 + 1)' 



(2) 



and the usual equilibrium cluster distribution / cq (n) 
given by 



/eq (n) = fo exp (-0G (n)) , 



(3) 



where G (n) is the work of cluster formation for a cluster 
with composition n = (m, n 2 ), we obtain 

df {ni,n 2 ,t) _ 
dt 

- k + (ni,n 2 ) / eq (ni,n 2 ) 

/ (ni,n 2 ,t) _ f{ni + l,n 2 - l,t) 

_/eq(ni,n 2 ) /eq (ni + 1, 7l 2 - 1) 

+ k + (m - 1, n 2 + 1) / C q (ni — l,n 2 + I) 
' / (m - l,n 2 + l,t) _ f{n 1 ,n 2 ,t) ' 
_/ cq {ni - l,n 2 + 1) f cq (ni,n 2 )_ 
a + (ni,n 2 )/eq (ni,n 2 ) 

f(m,n 2 ,t) /(m,n 2 + l,t) 



/eq(ni,n 2 ) /eq(ni,n 2 + l) 

+ a + (ni,n 2 - I) / eq (ni,n 2 - 1) 
f(m,n 2 -l,t) f{m,n 2 ,t) 



(4) 

/eq(ni,n 2 - 1) f cq {ni,n 2 ) 
which can be written in the form of continuum equation: 



dN (n) 
dt 



8 7 31 
dni dn 2 



-divJ = -V J, (5) 



where the components of the nucleation flux J are given 
by 



cq 1 9ni 9ti 2 f ' 



^n 2 



9ni 9n 2 



(6) 



(7) 
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and simplified notations 



(ni,n 2 ), a+ 



a + (n 1( n 2 ) and / eq = / eq (ni,n 2 ) and $ = $(n 1; n 2 ,t) 
are used, where 



$ (n 1 ,n 2 ,t) 



f (m,n 2 ,t) 
icq (ni,n 2 ) ' 



Apparently, 



<^fii ~l~ Jri2 ^ /cq 



dn 2 ' 



(8) 



(9) 



Therefore, the net flux coming into this composite nu- 
cleus is the incoming flux to the wetting layer of the in- 
termediate phase from surrounding parent phase whose 
magnitude is determined from the rate constant a + , since 
all materials are suppled from surrounding parent phase. 

Equations © and ([7]) can be put in the form of matrix 
equation using the short-hand notation d ni = d/dni and 
d 7l2 = d/dn 2 as 



which will be written in short 



-/ eq (n)R(n)V n <S>(n), 



(10) 



(11) 



where J n , <I> and V„ are row vectors, and R is a sym- 
metric square matrix defined through Eq. (1101) . The two 
flux J ni and J„ 2 are linked [18[ by the non-diagonal 
rate matrix R. By introducing the unit row vector 
= (e ni ,e„ 2 ), where the superscript T indicates the 
transpose vector, the nucleation current vector is written 
as 



J 



J - e J 



3n 2 Jn% 



(12) 



using the unit vectors e ni and e„ 2 along the Cartesian 
coordinate (ni,n 2 ). 

Since the Fokker-Planck equation given by Eqs. ©- 
([7]) for the growth of composite nucleus has the same 
form as that used to study the binary nucleation, we will 
extend the theory [27H29L [34M37I ] developed for the binary 
nucleation to study the scenario of nucleation and growth 
of a composite nucleus in the next section. 



III. CRITICAL NUCLEUS AND POST 
CRITICAL NUCLEUS 

A. Critical nucleus in the size and composition 
space 

In this section we extend the theory of nucleation flux 
of binary nucleation in the size and composition space 
developed by Fisenko and Wilemski [38[ to the composite 
nucleus shown in Fig. [T] for which the rate matrix R is 
non-diagonal. 



At this point we introduce a new coordinate, the size 
n and the composition x defined by 



n = m + n 2 , 
x = n 2 /n. 



(13) 
(14) 

Then the covariant formulation p39Ti4l| of the Fokker- 
Planck equation can be used, and the distribution func- 
tion is given by 



ip(n, x) = nf (n x ,n 2 ) 



(15) 



where the factor n comes from the Jacobian 
9 (ni, n 2 ) jd (n, x) = n. Then the Fokker-Planck 
equation ([5]) in the size-composition space is written as 



dip (n, x) 
di 



where the flux components are defined through 

J n = n(J ni + J n2 ) , (17) 

Jx = (1 ~~ ^) Jn.2 ) (1^) 

and Eq. ([lip is transformed into 



, Rx , n Rx , 



4>. 



(19) 



where the elements of the reaction rate matrix R defined 
through Eq. (fTS]) are given explicitly by 



Rn 



• ^ ' R%j — n Rtot, 



(20) 



Rn,* = Rx.n = -x (R n + R 21 ) + {l-x) (R 12 + R 2 {%) 

Rx,x = i (x 2 Rn + (1 - xf R 22 -2x{l- x) # 12 ) ,(22) 

using the elements of the matrix R, which will be written 
in short 



J X = -/eq (x) R {x)V x $ (x), 



(23) 



where J x and V x are the column vectors given explic- 
itly in Eq. flT9")) . Equation (|2"3")l is formally the same as 
Eq. (fTT]) . Therefore, various formulas for the steady-state 
flux in the binary nucleation can be used just by chang- 
ing the suffices n = (ni,n 2 ) to x = (n,x). Then, it 
is apparent from Eqs. (P2T)l) - (f2"2"|) that R n , n plays the role 
of Brownian diffusion coefficient in the size space, while 
R x ,x the Brownian diffusion coefficient in the composi- 
tion space |38| . 

Since R x , x — > as n — > oo, the fluctuation of the com- 
position x will be suppressed and x will become constant 
as the size of the nucleus n increases. The growth of nu- 
cleus is governed mainly by the reaction rate R n , n , and 
only the size of the nucleus increases. The coupling of the 
flux in size J n and that in the composition J x disappears 
when R n>x = 0, which is attained for the composition 



Xk 



i?12 + R 2 

Rtot 



(24) 
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where i?tot is defined through Eq. (|20|) . The diffu- 
sion coefficient in the composition space R x x given by 
Eq. ([22]) is also minimized for the composition Xk given 
by Eq. (|24l) . Therefore, the composition Xk is the kineti- 
cally optimum composition of nucleation, which is solely 
determined from the kinetic factor R. 

Next we will consider the nucleation flux at the sad- 
dle point (n*,x*) in the size and composition space and 
{n*,^) in the component space characterized by 



or 



dG 
dn 

dG 
dni 



dG 

dx 

dG 
dn 2 



0, 



0, 



(25) 



(26) 



where these two sets of equations are equivalent. 

From Eqs. (fTTj) and (|2"3")1 . it is obvious that we can use 
the results of the previous work [22j simply by replacing 
the matrix R by R. First, the angle cj for the direction 
of the gradient of $ in the size and composition space is 
given by the formula 



where 



tanw = s ± \/ s 2 + 1 



Rx,xG x x Rn.nGfi r> 



2 \ Rn,xG n ^ n -\- Rx.xGn 
Rn,nG n ^ x ~t~ Rn,xG x ,x 



Rn,xG n n 



R&.t.G 



x t x^ n,x 



and 



Gji : n 



d 2 G 

dn 2 ' 



Gx,x 



d 2 G 
dx 2 ' 



Gri x 



d 2 G 

dndx 



(27) 

(28) 
(29) 

(30) 



are the derivative (Hessian) of the free-energy G (n, x) 
in the size and composition space at the saddle point 
(n* , x*). These quantities are related to the derivative in 
the component space 



Gn — 



d 2 G 
dn 2 ' 



G22 — 



d 2 G 
dn 2 , ' 



G12 — 



d 2 G 
dn\dni 



through 

Gn^n 

G XiX 

Gn _x 



(1 - xf Gn + 2x (1 - x) G12 + x 2 G 2 2, 
n 2 (Gn - 2Gi 2 + G22) , 
-n (1 - x) Gn + n (1 - 2x) G i2 + nxG 2 
dG dG 



(31) 

(32) 
(33) 

2 

(34) 



dni dn2 

The last two terms in Eq. (|3~4"|) vanish at the critical point 
from Eq. (j26l . 

Equations (|2T1) - (|29I) are slightl y d ifferent in definition 
from those used by others [28|, [35[ since we have off- 
diagonal elements R n ,x = Rx,n- Among the sign ±, + 



sign must be chosen when R n ,xG n ,n + Rx,xG n ,x < oth- 
erwise — sign must be chosen. For a sufficiently large 
critical nucleus n* 3> 1, we have R n ,n ^> Rn,x 3> Rx,x- 
Also R n . x — when the composition is the kinetically 
optimum composition x* = Xk- Then, we have s 2 S> r in 
Eq. ([2T]). which results in to — 0. Therefore, the gradi- 
ent of $ (V;c$) will be nearly parallel to the size axis n 
when the size of the critical nucleus n* is large. There- 
fore, only the size rather than the composition is expected 
to increase. 

Using the result of the previous paper [22J], we can 
derive the formula for the angle </> of the direction of 
nucleation flux J x , which is given by 



tan <h = -=- '' 



R x x tan a; 



Rn,n + Rn,x tanw 



(35) 



Again, for a sufficiently large critical nucleus n* 1, we 
can approximate 



tan < 



Rn,x + Rx,x tana; 

Rn,n 



(36) 



and 4> — as R n , n 3> Rx,x- Then the nucleation flux J x 
also becomes parallel to the size axis n as expected. 
Finally, the nucleation rate is given by 



foe 



-f}G* 



\ 









Ai 


det.R- 






A 2 



(37) 



where Ai and A2 are the negative and the positive eigen- 
value of the matrix product GR where 



y"~f ( G n n G n $ 



(38) 



and R is defined through Eq. (|T9"|) . The matrices R and 
G must be calculated at the saddle point. Therefore the 
last two terms in Eq. (|34[) are zero. 

It can be shown directly by calculating the eigenvalues 
Ai and A2 that they are related to the eigenvalues Ai and 
A2 of the matrix products GR through 



A, 2 = nAi • 



(39) 



Also, it can be easily shown that deti? = detii. There- 
fore, the nucleation rate Eq. (13"T1) is also given by the 
original formula 



I = foe 



-PC 



IdetR 



Ai| 
A 2 



(40) 



in the original (ni,ri2) composition space as expected. 

Qualitative assessment of the relative magnitude of the 
composition fluctuation and the size fluctuation near the 
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saddle point can be possible. We follow the argument of 
Fisenko and Wilemski [38| and write Eq. ([TO]) as 



A$ - A$ 



CO I " "lb. lb A 

An 
A$ 



An 



A$ 



A$ 



«^ic — ./cq [ Rn.x a ~t~ Rx,x A 



— —fcqRx,x-£-^ (1 + 



An 
A$ 
Ax 



Ax 



where 



R x , x An' 



(41) 

(42) 

(43) 
(44) 



measure the relative magnitude of the cross terms in 
Eqs. (gB and (05). Since 



(45) 



An ~ v/l//3|5 2 $(n*,a;*) /<9n 2 |, 
Ax~ \/l//3|5 2 *(n*,x*)/9x 2 |, 



(46) 



are the curvatures of the free-energy landscape at the 
saddle point, and 



Rn,x — Rx,n — Riot \<^k ^) ) 



(47) 



we hnd 

Z: 

H '. 



(x k -x*) / d 2 $(n*,x*) /dx 2 



n* \i —d 2 <& (n* ,x*) / dn 2 ' 
n*R t pt (xk - x*) 

(1 - X*f i?22 + X 2 Ru - 2X* (1 - SB*) i?21 

l-d 2 <P(n*,x*) /dn 2 
V d 2 <I>(n*,x*) /dx 2 ' 



(48) 



(49) 



When Z <C 1 and H <C 1, the cross terms in Eqs. (|4"Tj) 
can be neglected, and we can study the time scale of the 
relaxation in size space r n by approximating Eqs. (|16[) 
and (JH]) by 



<9f t„ <9n cq <9n nAn 2 ' 



from Eq. (|45l) . which yields 



(51) 



for the relaxation time in the size space. Similarly, the 
time scale of the relaxation in the composition space t x 
is obtained by approximating Eqs. (|16p and (|41[) by 

dip A v -0$ i^,* 

ft ~ * fe /oqi? ^ ~ nA^' (52) 



which yields 



Therefore, the parameter 



H 



1~x Rn,nAx 

t„ R x j, An 2 Z ' 



(53) 



(54) 



will determine the relative magnitude of the relaxation 
times t„ and t x . If 



W 



(n*) i?tot 



(1 - .t*) 2 R 22 + (x*f Rn - 2x* (1 - X*) R 21 



-d 2 <5>{n\x*) /dn 2 
d 2 <$>{n*,x*) /dx 2 



(55) 



calculated from Eqs. (gSJ and (J49J) satisfies VF < 1, then 
*C t„, and the Brownian diffusion along the compo- 
sition axis is the fastest process and that along the size 
axis is the slowest process. Then the fluctuation of the 
composition will be suppressed and the nucleation flux 
will be narrow along the composition axis. On the other 
hand, when W 3> 1, the nucleation flux will spread along 
the composition axis even at the saddle point. In fact, 
since W 3> 1 in the model calculation by Wyslouzil and 
Wilemski IZfill which shows ridge crossing, Fisenko and 
Wilemski [38| argued that we can expect some effect re- 
lated to the ridge crossing when fF>l. 
From Eqs. (gT]) and we have 



J X 
J n 



(1 + H) Ax 
W{l + Z) An 



(56) 



and {Ax/ An) <C 1 in usual condition [38| when the sad- 
dle point is located within a narrow deep valley, we may 
expect J x / J n *C 1. Then the nucleation flux at the sad- 
dle point along the size axis is much larger than that 
along the composition axis. 

Returning to our original problem of composite nucleus 
for which the reaction matrix R is given by Eq. (I10[) , we 
find 



Rn,n — not , 
Rn,x = (1 - x)a + , 



1 



K + (1 — x) a 



2 „,+ 



(57) 
(58) 

(59) 



and i? t ot = o; + , and the kinetically optimum composition 
is given by Xk — 1 from Eq. (|2"4")l . We need information of 
the Hessian of the free-energy G to calculate the angle to 
of the gradient of $ from Eq. (|2"T1) and the angle 4> of the 
direction of the nucleation flux from Eq. (|35p . However, 
it is apparent that — > as n — > oo from Eqs. (|57l) - (|5"9")) 
and (|36p . Nucleation rate is proportional to I oc V a+ n + 
from Eq. (|4T)|) [22[. The time scale given by Eqs. (j5"Tj) and 
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(I53p can be written more explicitly by 



An 2 



k+ + (1-ar) a+ 



(60) 
(61) 



Therefore, the time scale of nucleation r„ in the size space 
is determined from the curvature An of the free energy 
surface along the size axis and the attachment rate a + 
from the parent phase to the intermediate metastable 
phase. The time scale t x in the composition space, on 
the other hand, is determined from the curvature nAx 
of the free energy surface along the composition axis and 
a + as well as the transformation rate n + from the inter- 
mediate metastable phase to the final stable phase. At 
the kinetically optimum composition x\, = 1, however, 
the size and the composition is decoupled and t x is de- 
termined solely from nAx and k + . 

Figure [5] shows an example of the free-energy surface 
calculated from our model [9j in the size n and the com- 
position x space. The size n is scaled by the typical size 
of the critical nucleus n* . There are two saddle points of 
almost the same size with different composition indicated 
by two points in Fig. [3J the one that corresponds to the 
critical nucleus of the metastable intermediate phase only 
(x = 1) and the other that corresponds to the core-shell 
nucleus (Fig.[T]) that consists of a core of the stable phase 
surrounded by a shell of the intermediate phase (x ~ 0.3). 
These two saddle points are not crossed sequentially by 
a single nucleation pathway. Rather these are crossed in 
parallel by two nucleation pathways shown by two heavy 
lines in Fig. [^independently. These nucleation pathways 
are the minimum-free-energy path (MFEP) [42| that cor- 
responds to the steepest-descent direction of the free- 
energy surface. 

In our model (U, figure [2] seems to suggest An/n ~ Ax 
since the size n is scaled by the typical size of critical 
nucleus n* ~ 10 — 1000, but we also need the magnitude 
of the attachment rate a + and the transformation rate 
k + to determine the time scale r n and t x of the nucleation 
quantitatively around these critical points from Eqs. (|6"0")) 
and (jSTj) . 



B. Post-critical nucleus in the component space 

So far, we have considered the critical nucleus at the 
saddle point. In order to study the nucleation pathway 
of not only the nucleation but also the growth of post- 
critical nucleus, it is more convenient to study Eq. ([FT]) 
in the original {n\,n 2 ) component space. This equation 
is rewritten in the form 



-f cq RV n $ = -RV n f + fh 



(62) 



where the first term in the right-hand side is termed the 
diffusion flux and the second term the drift flux. The 



0.8 



0.6 



0.4 



0.2 




0.5 



1 

n 



1.5 



FIG. 2. An example of the free-energy surface in the size (n) 
and composition (x) space calculated from the model based on 
the capillarity approximation The size n is scaled by the 
size of the critical nucleus n*. There exist two saddle points 
indicated by two points. The one on the x = 1 axis is the 
critical nucleus that consists only of the intermediate phase, 
and the other is the core-shell nucleus that consists of a core of 
the final stable phase surrounded by a shell of the intermediate 
phase shown schematically in Fig. [1] Two heavy lines are 
two independent nucleation pathways that correspond to the 
minimum-free-energy paths. 



latter is given explicitly by 

h = -f3RV n G, 



(63) 



which determine the growth of post-critical nucleus 
though this drift flux is negligible compared to the diffu- 
sion flux near the saddle point as dpi dm = 0. Eq. ([63]) 
is known as the Zeldovich relation [30, [3l[ . 
Explicitly, Eq. (|63p is written by 



(64) 
(65) 



However, as pointed out by Stauffer [35j, the continu- 
ous equation ([TT]) is originally derived from the discrete 
Master equation (j4|) of the attachment and the detach- 
ment of monomers, therefore we should use the difference 
equations 



dG 


r ti\2 q 

dn 2/ 


dni 


dG 


f R2 dn 2 ) 


dn\ 



dpG = J_dU 
dn\ 



df3G 
dn 2 



Icq 
1 - 
1 



dn\ 

dfeq 



Icq 

l - 



dn 2 

p P{G{ni,n 2 +l)-G{ni,n 2 )) 



(66) 



(67) 
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instead of the partial derivatives dG/drii and dG/dn 2 . 
Then, Eqs. (|64[) and (|65p are written more appropriately 
by 



(68) 



hi 


— Rn 


(- 


e P(G{m+l,n 2 ) 


-G(ni,n 2 )) 




+ R12 


(>- 


e j3(G(n u n 2 + l) 


-G(ni,n 2 )) 


h 2 


= R12 


(- 


e ^(G(m+l,n 2 ) 


-G(ni,n 2 )) 




+ R22 


fi- 




-G(ni,n 2 )) 



For sufhcicntly large clusters {rii — » cxd,?i 2 — > 00), we 
can neglect the surface tension of the spherical nucleus 
and the free energy of the nucleus is approximately given 
by [35[ G(ni,n 2 ) — —Afiini — Aa 2 n 2 , where A/ii > 
and Au 2 > are the chemical potential of the stable 
phase and the metastable intermediate phase relative to 
the metastable parent phase. Then, Eqs. ([68]) and (|69|) 
are approximately given by 



hi = R n (1 - e- pA ^) + R 12 (1 - e 



hi = R12 (l 



R22 (l - e 



0An 2 
/3A^ 2 



(70) 
(71) 



Therefore, the angle <f> of the direction of the nucleation 
flux is given by 



tan i 



«2 
"1 



-/3A(U 2 



^12 (1 - e-^ A ^) +^22 (1 

Rn (1 - e-^i) + i?i 2 (1 - e"^ A ^) 



(72) 



Since the number of molecule changes linearly in time t 
as rii = nit and n 2 = h 2 t, the hnal composition 



x 2 = 



n 2 



rn + n 2 



satisfies 



x 2 n 2 n 2 
= — = — = tan 1 

1 — X 2 Til Til 



(73) 



(74) 



given by Eq. (fT2"j). For large supersaturations A/ii — > 00 
and Afj, 2 — > 00, we find 



•1'2 



= tan ( 



I? 



12 



R 



22 



x 2 



R 



12 



Rn 



l-x k 



(75) 



Therefore, the final composition becomes the kinetically 
optimal composition x fc given by Eq. (|2"4"1) . 

Returning to our original problem with the reaction 
matrix R given by Eq. (|10p . we find Ri 2 + R 22 = a + 
and i?i2 + Rn = 0. Equation (l75t predicts that the final 
composition is x 2 = Xk — 1. Then, the post-critical nu- 
cleus consists only of the metastable intermediate phase 
and cannot reach the final stable phase with x = 0. 

These seemingly unphysical results are due to the as- 
sumption that both the stable phase and the metastable 
intermediate phase are comparably stable (A/xi — > 
00, Afi 2 — ¥ 00). In fact, we should have A/ii > A/12 ~ 0. 



Using the explicit form R in Eq. (|10[) , equations ((70)) and 
([7T|) are given by 



Til 



n 2 



-k+ (e-^ - e 



(76) 

~' 3A ^)+a+(l- e -' 3A ^)(77) 

Therefore, the number of molecules rii of the stable phase 
will increase as hi > 0, while n 2 of the intermediate 
metastable phase may increase or decrease depending not 
only on the chemical potentials A/zi and A/u 2 but also on 
the reaction rates a + and k + In the late stage of the post- 
critical nucleus when the depletion of monomer starts 
to occur, the incoming flux will stop and a + ~ 0, and 
the second term of Eq. (f77|) can be neglected. Then the 
number of molecules n 2 of the metastable intermediate 
phase will decrease to compensate the increase of the 
number of molecules ni of the final stable phase since 
h 2 = —hi < 0. Finally the supercritical nucleus will 
composed only of the molecules of the stable phase m, 
and the final stable phase with x — will be reached. 

In contrast to two minimum-free-energy paths indi- 
cated by two heavy curves in Fig. [21 the two real nucle- 
ation paths should gradually turn toward x = axis af- 
ter they overcome saddle points indicated by two points. 
They should merge into the x = axis as the size n 
increases according to Eqs. (|76]l and ([77]) since x = 
(n 2 — 0) in the stable bulk phase. This evolution of 
the supercritical nucleus will not involve saddle points 
and, therefore, must be a barrier-less gradual change. 
The time scale of evolution is determined mainly by the 
reaction rate k + and a + from Eqs. ([76)) and ([77|. Al- 
though both the free-energy landscape and the reaction 
rates play decisive role to determine the kinetics of nu- 
cleation at the saddle point, the details of the free energy 
landscape are irrelevant to the kinetics of the post critical 
nucleus. 

Of course, if the metastable phase is as stable as the 



final stable phase (A^i > Au 2 ), prediction of Eq. (f75)) 
suggests that the nucleus made of metastable phase (up- 
per nucleation route with x = 1 in Fig. [2]) could have 
a long life time since the kinetically optimum composi- 
tion is Xk = 1- Then the metastable intermediate phase 
would be macroscopically observable. The observability 
of the metastable phase around the core stable phase as 
the core-shell structure (lower nucleation route in Fig. [5]) 
also depends sensitively on the difference of the chemical 
potential A/ii and Aa 2 from Eq. (|77p . 



IV. CONCLUSION 

In this paper, we have studied the nucleation pathway 
of the critical and the post critical composite nucleus with 
core-shell structure not only in the size and composition 
space but also in the component space. By extending 
the results of the previous paper we could study 
the critical nucleus at the saddle point in the size and 
composition space. Our results suggested that the crit- 
ical nucleus can be more appropriately characterized in 
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the size and composition space. However, the kinetics 
of post-critical nucleus can be studied more easily in the 
original component space. 

Recent theoretical as well as experimental [H, 

UU results suggest that the appearance of the compos- 
ite core-shell nucleus does not necessarily mean two suc- 
cessive activations by crossing two saddle points sequen- 
tially [3, [ll|, [32J, |33[. Rather, the word two-step simply 
means that the nucleation pathway takes a roundabout 
course on the free energy landscape via the single saddle 
point which corresponds to the composite nucleus with 
core-shell structure. Our result in this paper will be use- 
ful to understand the nucleation pathway of such a two- 
step nucleation not only at the saddle point but also at 
the late stage of growth after crossing the saddle point. 

It must be noted, however, our analysis is completely 
confined to the steady-state process. The transient prop- 
erties can only be studied numerically by solving coupled 
Master Equations [l9j, |2jJ, |43| ■ In addition, our analysis 
is conducted under the assumption that the nucleation 
flux goes through the saddle point. Saddle point avoid- 
ance [13, l44U46| will be important if the anisotropy of 



the reaction matrix R is large or the ridge between sad- 
dle point is low, which can occur at high temperatures or 
near the spinodal point. In such a case, the ridge-crossing 
rather than the saddle-crossing may occur. Then the nu- 
cleation flux will spread over the whole phase space and 
the picture used in this study may break down. Also 
the growth of composite nucleus will be expected to be 
more complex. Finally, it should be noted that the recent 
progress to extract the reaction matrix R from simulation 
data [1^, [47J makes it possible to apply our formulation 
directly to a more realistic situation. 
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